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Abstract 



' In the context of Kaluza-Klein theories, we consider a model in which the 

^ , universe is filled with a perfect fluid described by a barotropic equation of 

r — ' state. An analysis of density perturbations employing the synchronous gauge 

o 



shows that there are cases where these perturbations have an exponential 
growth during a de Sitter phase evolution in the external space. 
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^ : 1 Introduction 

One of the main still open questions within the realm of modern cosmology is the 
origin of the primordial density fluctuations which played the role of seeds for struc- 
ture formation. Within the framework of gravitational instability, the two alternative 
scenarios which could lead to the generation of perturbations are either the amplifi- 
cation of quantum fluctuations during an inflationary era or fluctuations seeded 
by topological defects 0. Both inflation and topological defect models lead to ap- 
proximately scale invariant Harrison-Zel'dovich spectra on large angular scales, thus 
both are consistent with the COBE-DMR results P]. The characteristic features 
of the Doppler peaks in the angular power spectrum of the cosmic microwave back- 
ground, might be an important discriminating tool between inflationary fluctuations 
and those produced by topological defects H]. The dynamics of the initial density 
fluctuations has been studied by a large number of authors, either in the context of 
general relativity, or in multi-dimensional theories. 

In a generic inflationary model, the inflaton field (p experiences quantum fluctua- 
tions on all scales smaller than the effective particle horizon at the onset of inflation. 
During the era of exponential expansion, the scales inflate outside the horizon and 
therefore quantum fluctuations freeze in as classical fluctuations of 0, leading to gen- 
eration of energy density perturbations. The amplitude of the density perturbations 
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of a given scale A when they cross back the horizon at the end of inflation is [^] 




(1) 



Thus, through a linear mechanism one obtains a spectrum of density perturbations 
of constant amplitude, which under certain choices lead to the required amplitude 
of initial fluctuations, sufficient to generate structure formation without distorting 
the isotropy of the cosmic microwave background. 

For a number of authors, inflation is plugged with difficulties (e.g., potential 
fine tunning of coupling constants in order to achieve the required amplitude of 
density perturbations) and therefore we believe that it is important to investigate 
alternative scenarios. As such, one can study topological defect models (which will 
not be addressed here) or the growth of density perturbations in multi-dimensional 
theories, in particular, during a de Sitter phase. One may wonder whether in order to 
get non-zero density perturbations during a de Sitter phase while remaining within 
the framework of general relativity, one should either consider theories where the 
energy momentum is not conserved, or to remove the condition of a perfect fluid. 
We will sketch why none of these attemps work and therefore motivate our choice 
to work with multi-dimensional theories. 

A scalar-tensor theory of gravity which involved a non-conserved energy momen- 
tum tensor was proposed by Jordan but it was criticized [0, basically because 
it fails to successfully incorporate non-relativistic matter. A study of the evolution 
of scalar perturbations in a universe filled with a viscuous fluid, in the framework 
of general relativity, shows that the density perturbations decay (see, appendix A). 
So, to investigate whether scalar perturbations can grow during a de Sitter phase, 
we will consider a universe filled with a perfect fluid described by a barotropic equa- 
tion of state, in the context of Kaluza-Klein P] theories, where energy momentum 
conservation allows density perturbations provided the de Sitter phase characterizes 
only the four- dimensional subspace. 

In this work we are only dealing with the evolution of scalar perturbations, 
without examining their origin. In the context of inflation, for example, initial 
perturbations are generated by quantum fluctuations of the inflaton field. We will 
perform our analysis employing the synchronous gauge, while keeping in mind the 
residual coordinate freedom in order to finally keep only the physical modes. 

Our presentation is organized as follows: In section 2 we first set up the prob- 
lem, specify the notation, write down the background equations and the energy 
momentum tensor. In section 3 we introduce linear perturbations and derive the set 
of equations that describe the evolution of scalar perturbations in the synchronous 
gauge. In section 4 we solve the set of equations for two particular cases, where 
the external space is in a de Sitter phase. While our choices for these cases may be 
rather special, our aim here is only to show that it is possible to obtain growth of 
perturbations; the issue of how generic this situation is, is left for future work. We 
will show that, within the framework of multi-dimensional theories, we can indeed 
find at least one example where there is growth of scalar perturbations, in particular 
with an exponential behaviour, during a de Sitter phase. We summarize our results 
in section 5. We close with two appendices; in the first one we discuss density per- 
turbations in a viscous universe, while in the second one we analyze the problem of 
the residual coordinate freedom. 
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2 The model 



The Lagrangian density in a n-dimensional spacetime, coupling gravity to ordinary 
matter, is 

L = ^^^/^R-Lm, (2) 



where all quantities are defined in the multi-dimensional spacetime. From this 
expression, we can deduce the field equations in higher dimensions 



Rab — -QabR 



T 



AB 



StiGTab 
. 



(3) 
(4) 



We consider that each spatial section is divided into two spaces, the external one of 
dimension di and the internal one of dimension d2, such that n = 1 + cii + ^2- We 
assume the metric to be of the form 



-dt^ + a{tf'^ijdx'dx^ + h{tY'^abdx"-dx^ , 



(5) 



where is the external metric and 7aft the internal one, and both spatial sections 
have a constant curvature, not necessarily zero. Normalizing the curvature, we write 
fci = — 1, or 1 and ^2 = — 1, or 1 for the curvatures of the spaces with dimensions 
di and d2 respectively. 

We consider an anisotropic energy momentum tensor, with different pressures in 
each subspace. Thus, it takes the form 



T 



AB 



where 



with i 



,di', a = 1, . . . , (i2. 
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We consider that both pressures have a barotropic equation of state: 

Pi = "iP ; P2 = a2P ■ 



(6) 
(7) 

(8) 



The differential equations we obtain from Eqs. (|]) and 

a(t), b(t) and p are 



relating scale factors 



, a ,6 

"1- + "27 

a 



SttG 



di + d2 
d 



di(l + ai) + d2(l + 02) 



P 



- + (rfi-l)(-)^ + rf2-^ + (cii-l)^ 
a a a 



SttG 



di + d2- 1 



1 — ai + ^2(0:1 — 0:2) 



b , -,n/&n9 , dh _^.k2 

- + (rf2-l)(7)' + ^^l-7 + (4-l)T| 

ao 0^ 



(9) 



(10) 



SttG 



di + d2- I 



1 — a2 + di{a2 — ai] 



P 



d b 
p + di{l + ai)-p + d2{l + a2)-rp = 0. 
a b 

Some special solutions for these equations were found by Salidev |0. 



(11) 
(12) 



3 Perturbation equations 



We now proceed with the perturbative level. We introduce in Eqs. 
quantities 

'^gAB + hAB , p = °p + (5p , p = ^p + 5p 



and 



9ab 



) the 
(13) 



where ^Qab^ °P and represent the background solutions while Habi and 6p 
are small perturbations around them. We will also impose that the perturbations 
behave spatially like plane waves. Due to the anisotropy of the space, the equations 
take a tractable form only if the wave is defined in just one space. We consider that 
all perturbed functions depend only on the coordinate of the external space, i.e., 
S{xi,t) = 6(t) exp{iq- x), where Xj denote the coordinates in the space of dimension 
di. 

The equations are invariant under an infinitesimal coordinate transformation 



X 



^A _|_ ^A^ have the freedom of imposing a coordinate condition. We 



will impose the synchronous coordinate condition [1^ 

hAo = . 



(14) 



As it is well known, there is yet a residual coordinate freedom |11|, represented by 
a coordinate transformation that preserves the coordinate condition. This can lead 
to non-physical modes in the final solutions. We will discuss later this problem, 
establishing the criterium to eliminate them (see, appendix B). 
We define h = h^/a^., H = hj^/b"^, A = Sp/p. After a long but straightforward 
calculation, we obtain from Eqs. (^ - (0) the following system of differential equa- 
tions 



h + 2-h + H + 2jH 
a b 



H+{di- + 2d2j)H+ 
a b 



2do 



b / , , N /&n9 , db , , .k2 

- + f/2-1 7 ' + ^1-7+ ^2-1 tI 

ab ¥ . 



H 
A 



2(^1^ + ^27) A 
-d2Tn + 



(15) 



A + (1 + ai)6u\i + i(l + ai)h + i(l + a2)H 



[1 + ai)6u' + (1 + ai] 



d 6 

(2 - diai) a2Q;27 

a b 



6u' 



--(ai - a2)H'' - aiA-' 



(16) 
(17) 

(18) 



These equations describing the evolution of scalar perturbations, are too complicated 
to be solved in the general case. In the next section, we will solve them in two 
particular cases, where the equation of state in the external space is pi = —p. 
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4 Specific cases 

Here, we will solve the system of differential equations describing the evolution of 
density perturbations in the framework of multi-dimensional theories with a = —1, 
for the following particular cases: 

(i) Both the external and the internal spaces have scale factors with a power-law 
behaviour, a oc f (r > 0) and b oc f^, respectively. The internal space is also flat. 
To get a no n- vanishing 6p/ p, the set of equations describing the evolution of density 
perturbations requires 0^2 7^ —1. On the other hand, the background equations and 
the energy conservation equation imply 



^2(1 + «2) 

d2s{s - 1) 



+ di — 1 

One can easily check that Eqs. (p^Sf ) - (0) imply 

A = -^(1 + a2)H + const 
and lead to the following third order equation for H : 



(19) 
(20) 

(21) 



H'" + H" 



1 



ri(l — r) 



dir + d2S — r + 1 



+H' 



+H 



2 



{di{r^ + 2r) + d2s{2 + 3r - 2s) - - 2} 



+ 



ri{l — r) [ri(l — r)]^ 



{rfir^ + d2s{2r - s) - 2r + 1} 



0, (22) 



where primes denote partial derivatives w.r.t. conformal time rj, defined by adr] = dt. 
After some manipulations, we obtain that, as far as time evolution is concerned, 
Eq. (|22|) reduces to the second order equation 



d'^g 1 dg 
dr'^ T dr 



1 1 
724(1 -r)2 



{1 + r\di - 4di + A) 



+s'^d2{d2 + 8) - r(6rfi - 4) - 6d2S + 2rsd2{di - 6)} 



where 



1/(1-0 



r 



7 
2 



3 — dir — d2S 



.7 . 



2(r - 11 



2 2 
q 7] 



(23) 

(24) 
(25) 
(26) 



Equation ( [23| ) is a Bessel equation, whose solutions are Bessel functions of order 
d^s^ + 2dls^{s + 2) + djs^is + 4) + Sdjs^ + 4^2^(5 + 2) + 4 



A[d2s{2 -s) + 2]2 
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(27) 



This expression simphfies a lot if (ii = 3 and = —1, which means that the 
external space grows like with r > 1, while the internal one goes like where 
s < 0. In that case, 

-i±^ . (28) 



2(1 -r) 2 
The general solution of Eq. (p3[) is 

g{T) = uJ^{t) + vJ^J 



r) 



(29) 



where u and v are arbitrary constants. For the sign in the argument of the Bessel 
functions we choose r = +qr]. 



In terms of the density contrast, we find 
1 



A 



(1 + a2)ri-~^q^ j 7] (^uJ^{qr]) + vJ_^{qri)^dri. 



(30) 



For the particular case we are considering here, and re-expressing the solutions in 
terms of the cosmic time t, we find the asymptotic behaviour for A, for t — > 
{jj oo) and t oo (?7 — > 0), namely 







oo 



d2+2 



A ^ t (ci COS t '^■^+02 sin t '^■2), 



A 



-d2-5±2(d2+l) 
t 2d2 



(31) 
(32) 



So, initially the density contrast has an oscillatory behaviour with decreasing am- 
plitude, while asymptotically it tends only to decreasing modes. 



(ii) The external space is a de Sitter flat space filled with a perfect fluid described 
by Pi = —p, and its scale factor goes like a cx e^*, where r is a constant. The in- 
ternal space has dimensions d2 > 3, constant non-zero curvature and constant scale 
factor. The energy conservation equation implies p = const and the background 
equations lead to ^2 < 2/(1 + 02). To get a non-zero Sp/p, the perturbation equa- 
tions require 02 7^ —1. The system of equations describing the evolution of the 
scalar perturbations simplifies to : 



H + dirH - 



2{d2 - 1) 



h + 2rh + H 
H 



q_ 



A + -{l + a2)H 



Gr^A ; 

(rf2-l)^Arf2 



0; 
. 



Clearly, Eqs. (|| 



|36|) again imply 
1 



A 



[1 + a2)H + const 



(33) 
(34) 

(35) 
(36) 

(37) 



To get a growing solution of Eq. (^4]) we find the same requirement as the one 
imposed by the energy conservation equation, namely 



do < 



1 + a2 



(38) 
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Then, the solution for the density perturbation is 

A = f 2 ci J^(gf) + C2J-u{qf) 



(39) 



where f = 




and 




(40) 



The asymptotic behaviour is easily obtained. For small t, large f, 



A = e-5('^i-^)*cos(/3e-^* + 5) , 



(41) 



where /? and 6 are constants. For large t, small f, we obtain a mode which exhibts 
exponential growth, 



So, the density perturbations exhibt initially an oscillation with decreasing amplitude 
and then evolve towards a behaviour characterized by an exponential growth. 



One of the most interesting and important issues within modern cosmology, is un- 
doubtfuUy the origin of large scale structure. Within the framework of gravitational 
instability, there are two classes of theories attempting to answer the question of 
structure formation: inflationary models and topological defect scenarios. The cos- 
mic microwave background anisotropics, attempting to bridge between theoretical 
models and observational data, might support or rule out one of the two families of 
structure formation. In the context of inflation, quantum fluctuations of the scalar 
field give rise to density perturbations only once inflation ended and the perturba- 
tions re-entered the horizon. In this work, we are mainly interested in a model where 
density perturbations can grow during the inflationary era. We are working in the 
context of Kaluza-Klein theories, where the external space inflates. 

We employ some special solutions of multi-dimensional theories, and study more 
closely two particular cases, one with power-law behaviour of the scale factor, and one 
with exponential-type behaviour. The flrst one leads only to asymptotical decaying 
modes, while the second one, in spite of an oscillatory decreasing initial behaviour, 
gives asympotically an exponential growth of density perturbations. This indicates 
that such mult i- dimensional theories may furnish sufficiently large perturbations in 
the beginning of the radiation dominated era, thus they may provide an alternative 
picture. However, one should keep in mind that a further analysis, for example 
the determination of the amplitude of quantum fluctuations, which may allow a 
confrontation with observational data, is required before a flnal statement is reached. 



Appendices 



A Perturbations in a viscous universe 



Another way of obtaining a growth of density perturbations during an inflationary 
phase, is to use an imperfect fluid in the energy-momentum tensor. This can be 



A 




(42) 



5 Conclusions 
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achieved, for example, by introducing viscosity. Here we consider bulk viscosity 
in the context of the four- dimensional general relativity. The bulk viscosity effects 
can be deduced just by replacing p i— >■ p — x(p)0 in the energy-momentum tensor. 



Eq. (I), where 6 
The field equations are 



We consider a barotropic equation of state, p = ap. 



3{-f 

Qi 

a a 



8nGp ; 



-8-KGap + 2ATcGxip)- ; 

a 



p + 3(1 + a)p = 9xi-y 



(43) 
(44) 
(45) 



If x(p) = Xo = const, we have the simple solution a oc e^*, where H = {1 + a)p/3Xo 
and p = const. If on the other hand, x(p) = XoP^ have the solutions found by 
Murphy [|l^], which reduce either to the previous one in the limit t ^ — oo, or to 
the flat perfect fluid solution in the limit t — oo. 

If we perturb Einstein's equations in the presence of viscosity, we flnd the fol- 
lowing coupled differential equations 



h + 2-h - UnGxih - ^) = ; 
a 

A + 9(-)2(^ - x')A + (! + «- 6^-)(^ -h) = 0; 



l + a-3--^- 
a p 



a p 

;i+«)(2 



p a 



a 

3a)- 
a 



a 



3{- + {-y +9a{-y- + 
,a a / a p 



a \ a p 



^ _ 



,a 



aA - 3x'-A 
a 



P 



, 



(46) 
(47) 

(48) 



where A = 6p/p, = 6\i and x' = dx/dp. 

We will now determine the perturbed solutions for the density contrast. If x = x^, 
x' = 0, then using the background relation 



l + a = 3{x/p){d/a) 



(49) 



we obtain 

A = . (50) 

This solution cannot be eliminated by a residual coordinate transformation. So, if 
the viscosity coefficient is constant, we have a de Sitter phase, during which the 
density perturbations do not vanish; however they decrease exponentially with time. 
Let us now consider the case x = XoP- t^e asymptotic limit t — >• —oo, the 
density contrast decays exponentially as the inverse of the cubic power of the scale 
factor at that stage. We find again a de Sitter phase, with non-vanishing density 
perturbations, which however decrease exponentially with time. On the other hand, 
in the limit t ^ oo, the perturbed solutions approach the ones found in the case of 
a perfect fiuid, like the background solutions do. 
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B Residual gauge freedom 



We consider the infinitesimal coordinate transformation x^ + x^-i under which 

the components of the metric tensor transform to 



QAB = QAB + X{A;B) 



(51) 



Imposing that the synchronous coordinate condition be preserved, we obtain the 
following solutions for the time- and space-components (in each of the two subspaces) 
of : 



X' = *(x)'^/^ + C(x)'\ 
r df 

x" = n^rj^+E{xr, 



(52) 
(53) 

(54) 



where C and S are arbitrary functions. Since the traces of the perturbed metric 
tensor, in each subspace, are 



2m(x) 



dt 








-2d,^^! 




a . 


dt 






fe2 




- 2(^21* 

b . 



(55) 
(56) 



one can show, using the definitions h — hj^ /a^, H — h^/h"^ and the equations for the 
perturbations, that the density contrast transforms under this residual coordinate 
transformation as 



A = A 



(1 + ai)di- 



a2)d2 



(57) 



Thus, any perturbed solution with a time behaviour {l + ai)di{d/a) + {l + a2)d2{b/b) 
can in principle be eliminated by a coordinate transformation and has no physical 
meaning. 
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